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Abstract 

An optimal effective current operator for flux qubit has been investigated by taking account of 
the inductive effects of the circuit loop. The whole system is treated as two interacting subsystems: 
one is the inductance- free flux qubit consisting of three Josephson junctions and the other a high 
frequency LC-oscillator. As the composite system hardly affords one excessively high energy LC 
photon, an effective theory for the inductive flux qubit providing its physical variable operators has 
been achieved, which can take account of the inductive effects but does not include the additional 
degree of freedom for the LC-oscillator. Considering the trade-off between simplicity and accuracy, 
it has been revealed that the optimal effective current operator resulting in an error only on the 
order of L^/^ provides an approximation of high accuracy, which is also verified numerically. 

PACS numbers: 03.67.Lx, 85.25.Cp 



1 



I. INTRODUCTION 



Superconducting circuits are promising candidates for quantum information 
processingii^i^ and, in order to reduce the impact of both charge and flux noise, flux 
qubit consisting of a superconducting loop interrupted by three Josephson junctions 
(3jj) has been proposed, designed and realized.-i^'^i^'^i2i»^'^ The loop in the original 
design is small enough and its inductive effects, therefore, could be neglected at the first 
approximation.- The constraint of the flux quantization on the three phases across the 3jj 
yields two independent phase variables for the system. On the other hand, inductive effects 
are essential in several inductive coupling schemea^^^^'^^^^^'^^^. These systems can be 
systematically studied by applying a general network graph theory^^ii^, which indicates that 
an independent phase is associated to the loop self-inductance in the original circuit and 
the 3jj flux qubit, thus, turns out as a three-phase system. In order to include the inductive 
effects judiciously, appropriate terms could be reallocated to improve the original operators 
in the two-phase system. First, the inductive effects, considered as corrections to the energy 
levels of the two-phase system, have been addressed but with some flaws by Crankshaw 
et. al. in a semi-classical approacb^S^ and, consequently, an effective Hamiltonian has been 
reached^ as well as a current operator in the two-flux-state basis for the flux qubit^. 
Another reason why we build up an effective theory to include the inductive effects is that 
an inductance of a non-negligible size may lead the device to a less useful qubit. 

Current operator is crucial to the accurate control, coupling and measurement of flux 
qubits.— In particular, it could play a key role in understanding the dynamics of the flux 
qubit by a general multilevel modeL-i^i^^ Although various forms of current operators have 
been utilized in all kinds of regimes, the validity of the specific current operators has not 
been justified seriously and error analyses are hardly available. In this work, a systematic 
investigation on the optimal two-phase effective current operator for the three-phase system 
is carried out and an error analysis is provided. 

The paper is organized as follows. In SecHIl we review some basic ideas on the loop 
current in a classical circuit model. In Sec lIIIt we construct the three-phase Hamiltonian 
and decompose it into a form showing that two subsystems weakly interact with each other; 
then we develop an effective theory, the photon transition path (FTP) approach based on 
the Brillouin-Wigner expansion^^, to describe the three-phase system in Sec JIVI In SecJV] 
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FIG. 1: (a)circuit of an inductive flux qubit witli tlie pliase difference (j) across the loop inductance 
L, the reduced applied external flux (px = 27r<I>x/*^*o with <I>o the flux quantum, and the phase 
difference (j)k across the A;th junction characterized via the critical current Ick and the capacitance 
Ck for k=l,2 and 3; (b) transformation between the current and voltage sources, the arrow and 
the plus/minus symbols indicate the directions of the current and voltage sources, respectively. 

we obtain the optimal effective two-phase loop current operator from the unique one for the 
whole system, and a brief numerical discussion is presented in SecJVIl 

II. CLASSICAL ANALYSIS 

The schematic circuit for the 3jj flux qubit with a loop inductance is demonstrated in 
Fig. [U^a), where the 3rd junction is a little smaller than those two others; representing the 
relative sizes, the parameter ak as 

= ^ = ^,A; = 1,2,3 (1) 

-'CO (-^0 

indicates the area factor of the kth junction, where Icq = (-^ci + ^02) /2 and Cq = 
(Ci + C2) /2 are design parameters. Parameters ai and 02 are supposed to be close to 
1, the deviations of which are determined by the accuracy of fabrication, while 0:3 ~ 0.8 and 
the reduced applied external flux (px is biased on the vicinity of (px = tt, all of which are 
selected to benefit the energy levels of the flux qubit. These three junction phases 0i,2,3 and 
the phase difference across the loop inductance L are not independent of one another and 
obey the flux quantization in this superconducting system as 

01 + 02 + 03 = + 0X, (2) 
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the signs of which are also indicated in Fig. [T]^a). 

In the classical regime, the junction performs as a current-fiux 2-port circuit element, 
a nonlinear inductance, and the flux quantization condition imposes a predetermined con- 
straint. In the DC regime, without considering capacitances, the loop current flows equiva- 
lently through four current elements in the loop including the Josephson junctions and the 
loop inductance as 

= /cfcSin0fc = = 1,2,3, (3) 

where 0i, 02, 4>3 and cj) are the possible static phase values obtained from Eqs. ([2]) and ([3]). 
Two opposite current directions present an additional degeneracy of the circuit. Further- 
more, in the AC regime, if only taking account of the small oscillations in the circuit, each 
junction works at the static phase point as a pure inductance 



$0 



2ttIc^ cos (j)k 

if cos(f)k 7^ for k=l,2 and 3. The series impedance of the circuit 

Q 

-iuL 



(4) 



Z(uj) = —iuo 



— iUUlJl. , ^ 



k=l 

with i the imaginary unit, provides its several characteristic frequencies; especially, when the 
circuit works at such an ultra-high frequency that the junction inductances can be treated 
as open circuits, there exists only one significant oscillation along the loop between its small 
inductance L and series capacitance Cger with a high characteristic frequency 



where Cser = (ELi C'fe ^)"^ = «serCo. 

Generally, the nonlinear effects of the junctions generate current components of new 
frequencies different from the external flux-driven source's. Only considering the output 
profiles of the junctions, we can still apply this kind of speciflc current sources to the rest 
of the circuit which obeys the linear superposition rules. Picking an arbitrary frequency lj 
in the frequency domain and utilizing a source transformation shown in Fig. [I](b), we have 



where Ick{'^) is obtained from Ick sin0fc via the Fourier transform. Interestingly, when L is 
small enough to neglect, Iioopil^) i'^ Eq. ([7]) does not depend on oj explicitly and we utilize 
the inverse Fourier transform JF~^ again as 

hoop{t) = ^ (-^«oop('^)U2L^o) = C'ser -^^^^7^—^, (8) 

k=l 

where 4>k = 4>x since vanishes when L ^ 0. This form of the loop current hoopif) 

directly goes with the fact that the junctions connect to a topological network consisting 
of linear circuit elements. Delightfully, hoopif) in Eq. ^ is in exact agreement with the 
one for the two-phase system derived in the quantum regime by Maassen van den Brink — 
and with our following effective one. This suggests that quantum superconducting circuit 
analysis and design might benefit in elegant ways from classical circuit theories and CAD 
tools. 



III. QUANTUM ANALYSIS FOR SYSTEM HAMILTONIAN 

To construct the Hamiltonian comfortably, we firstly select three junction phases 01,2,3 as 
the spatial variables and express the system Hamiltonian in a sum of energy terms similar 
to other superconducting loop circuits such as the RF-qubit and the SQUID-qubit as 




3p 



where Qk is the charge operator conjugated with the phase 0^, i.e., Qk = 

<pk, Qk = 2ei with e the electronic charge; Ej^ = Josephson energy of the fcth 

junction. The former sum in Eq.Q represents the total energy of junctions including their 
charge and Josephson energy and the latter term the loop inductive energy. According to 
the design, the reduced inductance size /3, 

= (10) 

is usually small enough that the loop phase difference behaves as a small variable with its 
norm ||0|| tending to be equal to zero, while the loop current still keeps finite due to the biased 



junctions. Consequentially, its conjugate variable Q^, which we refer to as 0, Q 



2ei, 



diverges on its norm according to the Heisenberg uncertainty principle ||0|| • ||(5</-|| ^ 



In the classical regime, a quadratic potential (^) |^ means that there is a generahzed 
restoring force 

^°v,,,„,30^ = M[i,i,ir (11) 



AttL ^ 2ttL 

providing a non-parallel generalized acceleration 

5o. = C-/„„=^f, (12) 

where Cdiag is a diagonal matrix with its diagonal elements being Ci,C2 and C3, and 

Oil «2 tts 

In the quantum regime, the deep quadratic potential explicitly in proportion to 1/(3 is 
capable to bind up the quantum states of this three dimensional system in the vicinity of 
a phase plane = 0, where a fast vacuum fluctuation occurs along the unique direction 
f parallel to the acceleration Qosc- Therefore, the original spatial variable set (</>!, 02, ^s), 
although helpful in the construction of the Hamiltonian, presents difficulties in handling the 
charge operator Q^^ which represents one of the most important quantum properties of the 
three-phase system. 

To solve the problem, we utilize a linear transformation to achieve another set of coordi- 
nates (0,^1,^2) where besides the other two coordinates are labeled via O-i and 82 and their 
conjugates are Qe^ and Qe^ , respectively. The linear transformation between these two sets 
of coordinates is introduced via a matrix A defined as [9i, 62, 4> + 4'xY — 02, ^s]"^, or 

equivalently as 

Qg = Qei-,Qe2-:Q4> = Qi:Q2,Q3 A ^. (13) 
Thus, the Hamiltonian H^p changes to 

Ha = ^QSAC^-^A^Qo + V{e,, 02, 0), (14) 

where V{6i, 62, 4>) is the potential in the new framework. Since the charge operator Q^f, tends 
to diverge when P —>■ 0, ii the charge coupling coefficients in AC^^^^A^ are assumed to be 
finite and independent of /3, a proper candidate for ©-subsystem on (^1, ^2) should avoid any 
direct charge coupling from the 0-subsystem. It mathematically requires that the directions 
of 9i and 62 in the original coordinates should be perpendicular to the acceleration direction 
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r of the oscillation mentioned above, which means that the plane spanned by 9i and 62 is 
unique as well as 

3 



Q4 



C ser ^ ^ 



fc=l 



Qk 

Ck ' 



(15) 



revealing the charge in the series capacitor C^er- Some other explanations in the classi 



cal regime are also given in Refs. 



14 



and 



201 . both of which have achieved the proper 



variable transformations by avoiding the cross charge energy terms between B and 
subsystems.— Although they have predicted the right ones based on the linearity of the 
circuit, it is more comfortable in the quantum regime to emphasize the reason why the G- 
subsystem as well as should be selected uniquely, since the diverging charge fluctuations 
merely serve as a pure quantum phenomenon. 

The remaining degrees of freedom endowed by A involve the internal variable selections 
of G-subsystem. A straightforward way is that 9i and 62 only deviate slightly from 0i and 
02, respectively; then, the whole transformation reads as follows. 
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-1 



(16) 



where the last term of its right-side is a set of constant biases as a translation in the 
superconducting phase space. For short, it can also be reformatted as 



(17) 



which clearly shows that 0^ reduces to 6k when 0^0. The transformed charge operators 



(18) 



indicate that Qe^ states the charge of the island between the junctions 1 and 3, and analo- 
gously for Qg^. If we also define 
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^3 



(19) 
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equating to (px — 0i — 02, three new phase variables ^1,2,3 confined by the fiux quantization 
seem to act as the junction phases in the two-phase system, which is confirmed by the 
following transformed Hamiltonian 



tr 



Ho + {a'^a + -] hwLc + Hint, 



(20) 



where 



1 ^ 
i^o = -Q^C^-'Q^^-J^^jfcCos^, 



k=l 



a) a 




1 

2' 



int 



3 3 / 

Ejk COS 6k Ejk COS i 6k + 

k=l k=l ^ 



Co, 



2p 



Ci + C3 C3 

C3 C*2 + C3 



(21) 
(22) 
(23) 

(24) 
(25) 



The Hamiltonian has been decomposed into three parts. The first part Hq is the Hamiltonian 
of the 6-subsystem for the inductance-free fiux qubit.- The middle part (d"'"a + l/2)hu!Lc 
shows that the LC oscillator consists of the phase variable and its conjugate and its 
characteristic frequency ujlc=^/ VLCser is high enough as mentioned in the classical regime. 
The operators and a are defined respectively as the photon creation and annihilation 

.2TT J Lh? ("^+«) 



operators and the dimensional factor of 



-^^Y c ~72~ proportion to P^^"^. The 

last part Hint is the interaction Hamiltonian between these two subsystems, which is weak 
enough compared with huJico^ci and Hq to make the G-subsystem only be slightly perturbed 
by the LC oscillator. The current implementation of Hint is useful for numerical solutions 
of diagonalizing Kronecker product matrices with the FFT tools,— which are also utilized 
in this paper, and its series expansion on 



H 



int 



k>l 



(26) 



where = ■^^-^mt|(/)=o, is fit for the perturbation methods which should deal with the 
couplings of different strengths. 
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Qubit 



\M2 



LC Oscillator Whole System 

FIG. 2: Energy diagram of flux qubit witli a loop inductance. When the inductance-free flux qubit 
and the LC oscillator interact with each other in a perturbation condition, the lowest eigenstates 
in the dressed state manifold Mq denoted with the dashed-line box are well separated from the 
ones in other manifolds TWi, • • • due to the large shifting caused by the LC-photon energy 



n 



FIG. 3: (Color online) Two typical photon transition paths (1) and (2) represented by the linked 
operator chains Hs,s{s)Hsfl and Hs,s{£)'Hs,sJ~Lsi,sA^)'Hs^fi, respectively. 

IV. PHOTON TRANSITION PATH METHOD 
A. Dressed states in manifolds 



To understand the energy diagram of this system shown in Fig. [21 let us briefly recapitu- 
late the well-known dressed state concept^^. For the sake of simplicity, we do not explicitly 
consider the possible inner degeneracy in the two-phase subsystem further and have a set of 
complete orthogonal basis Bfj^ = | Ho\^Pf) = ^f\Vf) | , where the arbitrary normalized 
eigenstate ly?/) goes with its eigenenergy Sj. The LC-oscillator keeps its sth eigenstate \VLs) 
as a)a\VLf) = s\Qs). 

When the interaction H^nt is neglected at the first approximation, it is convenient to find 
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that a series of dressed artificial-atom states {Iv^/) ® \ fio), ■ ■ ■ , \ipf)^\ftN), ■ ■ ■}, where 
is a non-negative integer, are the eigenstates of the whole system with their eigenen- 
ergy being {^f + \f\hJLCi ■ ■ i Sf + {N + ^)hujLc, ■ ■ ■ }, respectively. Since Hulc 
is much larger than ej independent of (3, the dressed states are so well separated from 
one another that the tensor-product states which keep the same photon number can 
be grouped together to form one so-called manifold. For example, the A^th manifold 
•^N — {Iv) \\Vf) ® I^Af))^lv'/) ^ ^Hq} consists of all possible eigenstates possessing N 
LC photons and maintains the same energy level structure as the two-phase flux qubit's if 
(A^ + is subtracted. After the weak interaction Hint turns on in the order analy- 

sis, the possible intra- and inter-manifold photon-assisted transitions bring perturbations of 
different strengths, which cannot completely destroy the manifold structures, so the per- 
turbed eigenstates in the A^th manifold n can still be distinguished from other manifolds' 
due to the (A^ + ^)hLL!Lc energy shifting. This kind of understanding can be revealed by 
one well-known perturbation approach, the unitary transformation (UT) method:^^, which 
introduces a specific unitary transformation T = e*'^|^^^t; resembling a time-evolution op- 
erator, to rotate the Hamiltonian Htr into a new one = T^HfrT so that it can be 
diagonalized as = Y1'n=o^n\^n){'^n\ on an arbitrary order of (3. The two-phase 
Hamiltonian Hn performs as an effective one for the A^th manifold AiN'- with the eigen- 
state basis Bj^^ = ||<^Ar) HnI'Pn) = ^nI'^'n) |, the A^th manifold A4n can be rewritten as 
Aijy = \(p) = T {\(piy) (g) \Qn)) ,^\'Pn) ^ ^Hj^^ ■ particular, when Hint is neglected 

we can select S = and obtain the effective Hamiltonian H^^ ~ (| + N)^lc + -^o indi- 
cating that the manifold M.n unsurprisingly becomes M.f^ when the interactions turn off. 
With the unitary operator T, one can also consequently construct other effective operators. 

On the other hand, since there is hardly an experimental way to keep the high-energy LC- 
oscillator excited in the superconducting circuit applications, what needs to be focused on 
actually is the lowest eigenstates belonging to the manifold M.^. This physical requirement 
also enables us to circumvent the additional discussions on that the inductance-free flux 
qubit as an infinite-level system still leaves the high-energy eigenstates in A^o not being well 
separated from but overlapping with the lowest ones in A^i for a specific value (3 in the 
energy diagram. By means of the Rayleigh-Schrodinger (RS) expansion with the arbitrary 
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eigenstate \ip) in A^o and its eigenenergy e being respectively expanded as 



(1) 



r(2) 



(27) 
(28) 



where both \^^'^^) and e*-'^'-' are in proportion to jS'^^^ for A; as an integer and \^^^^) belongs 
to the manifold M.'^q \ an effective Hamiltonian on the order of j3 has been obtained but 
without further discussions on the higher order expansions in Ref. [ij. In this paper, based 
on the Brillouin-Wigner(BW) expansion, another famous perturbation theory, we develop a 
photon transition path method to further explore the perturbation procedure and compare 
it with the one in Ref. 14 and also with the UT method. 



B. Formal substitution derivation 

Besides the order expansion in Eq. (1271) . we also expand \(p) in the energy eigenbasis of 
the oscillator as 



(29) 



s=0 s=0 A:=0 



where \ips) = (^s|v^), and {(pi ) oc /J*'/^. Since the expansion in Eq. begins with a 

constant number ^Hulc, we subtract it from Htr and redefine the Hamiltonian Htr as 



Hi 



tr 



Ho + hwLcCi^O' + H, 



'■int • 



(30) 



Consequently, the biased eigenenergy e satisfies e <^ fv^LC- The Hamiltonian Htr can be 
expanded clS cl Kronecker product matrix 



Hu- 



Hofi Hqi 

Hin -^1.1 + hiOLC 



H. 



HsA 



Ho,s 
Hi^s 

: Hs,s + shujLC 



(31) 



where the operator Hs^s = Hq + {Qs\Hint\^s) refers to the self-transition of the sth level and 
Hs,si = {^s\Hint\^si) = Hl^ g the transition between the sth and sith levels of the oscillator. 
The operator Hg^s consists of terms of different strengths due to the nonlinearities included 
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in Hinti and its dominant term Hq approximating Hg s on 0{(3^) suffers from an error on 
0(/3^/^) instead of 0(/5^/^) thanks to the optical selection rules. For the same reason, the 
dominant terms of the operators Hs,si are also weakened on ©(/J'*"*^'/^). More details about 
the order discussion are presented in Appendixes |X] and [Bl The eigen-equation Htr\^p) = e\(p) 
is decomposed into a series of equations as 



(^0,0-^) I'/'o) = -J^i^o.fcl^fc), (32) 
[shwLc + Hs,s - \ips) = - ^ Hs,k\'-Pk) (s > 0). (33) 

k^s 

The shifting of ^Jtwlc, thus, distinguishes Eq. fl32l) from the others in Eq. fl33|) . It is clear 
that when the loop inductive effects are totally neglected, this equation is capable to degrade 
into a two-phase eigen-problem as 

HoM = ^M), (34) 

and all of the other projected states \y^s^^){s>o) equate to zero, which suggests that the 9- 
subsystem decoupled from the LC-oscillator becomes an inductance-free two-phase system 
and there is no LC-photon excited at the first approximation. 

Since the LC-oscillator is of high energy, sJtlolc always dominates in Eq. (l33!l at the excited 
levels (we assume that the integer s is larger than zero in the following sections). To figure 
out the relative strength, a new set of operators are defined as 

Hs,. = (35) 



with an introduced Kronecker delta function 6a,b, thus Eq. (13311 yielding 

(l - HsA^)) \Vs) = Y^UsAVk). (36) 

k^s 

As the dominant term of 'HsAA is on 0(/3^/^), we expand the modified BW resolvent 
operator (1 — UsAAA^ ^.s X]j>o '^s s(^)' ^^"^ state \'^s) is given as 

i>0 siT^s 

This equation indicates that Iv^s), different from |(y9o), is a result of transitions from all of the 
other levels. Moreover, it is found that the large photon energy sItwlc and the corresponding 
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multi-photon processes involved punish all of this kind of transitions via diminishing them 
on specific orders of jS. 

According to the difference of the state \(Pq) from the others, we leave the right hand side 
of Eq. fl?r|) separated as 

\^s) = Y.Ksi^)^s,o\^o) + Yl E ka^)Am, (38) 

i>0 i>0 si^{s,0} 

which involves two types of PTPs to the sth level: the operator Ti* ^{e)7isfi means that the 
state \!fo) transfers from the ground level, then through arbitrary times of self-transitions, 
to the sth level and HI g{e)'Hs,si refers to the other state \ipsi) (neither \(po) nor l^^^)) from 
the Sith level. Since Si 7^ in the above sum, we also have 

n>o ii>Os2^{si,o} 
where s, Si and i have been substituted by Si, $2 and zi, respectively. The latter type of 
PTPs in Eq.f l38l) . therefore, can also be divided again as 

Ws) = Y.'^ls{e)nsflW^) + E Ksi^)^s,sM,,s,{e)ns,fiW^) + (40) 

i,ii si^{0,s} S2^{0,si} 

These substitutions employ the procedures of the BW perturbation approach with a clearer 
view on the orders of the terms on j3. For example, without the expansions of the BW 
resolvent, Eq. fHU]) resembles a familiar BW expansion as 

\^s) = 7— ^s.olv'o) + (41) 



' s,s 



1 - 1 - 

E 1 ; — ^si,ol</'o) + 

E E 1 ^Hs,s^ 7 ^ Hs,,S2\'^S2)- 

Furthermore, we are also able to substitute the 3rd part in Eq. fHUl) and divide it into two 
parts, the latter one of which can be substituted again. After retaining the transition paths 
to the sth level from the ground level and continuing this kind of substitutions for n — 2 
times with S3, S4, ■ ■ ■ , and s„ being introduced, we transform Eq. fl38|) into 

|^.) = Gi")(e)|^o)+E^il(^)lv^-)' (42) 
13 



where G^r\e) refers to all of the PTPs from the ground level involving no more 
than n times of non-self transitions and G^J^]„{e) all of the PTPs in a form like 
HI g{e)Ti.s,si'Hl\ si{s)i-isi,s2 ■ ■ ■ '^s"!! s„_i(^)'^sn-i,sn) ^he dominant terms of which are at least 
on 0{(3^^'^) contributed by Hujlc- Since n can increase so large as to make G^J^s„{£) negligible 
on an arbitrary order of f3, the sth level is uniquely determined by the state {(fo) with a 
corresponding operator Qsi^) being defined by 

\^s) = gs{e)\^o), (43) 

where Qo can also be added as an identity operator with l2p\'^o) = Iv^o)- For example, 
the operators ^1,2,3 (f^) are approximately given in Appendix [Bl Therefore, the projected 
state lifo) with a map, which a series of operators Qs{^) function as, covers the three-phase 
state \(f) completely and accurately. Equation (H3|l mathematically describes one physical 
understanding that for the states in the manifold A4q photons persisting in all of the excited 
levels come from the ground level via all possible PTPs as illustrated in Fig. [3] due to the 
perturbations of the inductance-free flux-qubit. 

Substituting \lPs) in Eq. (1321) with the aid of Eq. ( l43l) . we have an eigen-like problem 

Hie)\ipo)=e\ipo), (44) 

where the 2)seu(/o-Haniiltonian H{e) is defined as 

H{e) = Ho,o + Y,^oJs{e). (45) 

s>0 

In the definition of H{e), all of the terms in the latter sum can be described in a general form 
Ho^s'Hl^si^)'^s,siHl\^sii^)'^si,s2 ■ ■ ■'^s"Ii,s„_i(^)'^sn-i,o wWch can be interpreted in the terms 
of the photon-assisted transitions as that the LC photons spread to one specific excited level 
such as the sth one from the ground level through an arbitrary PTP (the role Qsi^) plays) and 
then return back ( an operator Hq s closes the whole PTP). Therefore, the PTPs introduced 
by the operator H{e) are not only linked but also closed, starting from and ending with 
the ground level. It should be emphasized that a one-to-one correspondence is established 
between the terms in this sum and the closed photon transition paths(CPTP). Putting aside 
the details of the CPTPs in this section, one idea can be accepted that the longer path the 
photons travel along, the weaker effects are brought. Based on the BW expansion, the above 
derivations do not lose any accuracy thanks to the formal substitutions we utilize. Yet, as 
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drawbacks, to make the whole solution available, we still need to deal with the infinite terms 
included in H{e) and its dependence on the eigenenergy e which is actually unknown before 
we successfully solve the problem. 

One common solution for these two problems is to employ the standard RS perturbation 
method, which utilizes the expansions of e and |<y9o) in Eqs. (1281) and (12^ . respectively, 
and all possible results can be achieved by checking terms on the same order of f3 in Eqs. 
fH3|) and (jS]). This approach mixes up the BW and RS perturbation methods and benefits 
at least on two aspects due to a fact that perturbation effects of different strengths are 
able to coexist in one photon-transition matrix element which we can manipulate in a more 
physical manner. One is that instead of the step-by-step style we directly expand Eq. 
on a specific order of P^^^ and, consequently, achieve a series of equations including all of 
the cases below this order. In this context, our method now acts as an improved wrapper 
for the order analysis utilized by Ref. [ij, and the difference is that we use the projections 
before the order comparisons while they prefer that the latter one goes first. The other is 
the convenience that we can more easily predict characteristics of the perturbation results. 
For example, without the emphasis on sfuvic and the consequential result Eq. ( H3i) . it is not 
obvious in the previous paper that the projected states on the excited levels can be derived 
from Iv^o), although the term sHujlcI'Ps''^'^'^) in the expansion of Eq. ([33D on 0(/3'=/^) with 
k being an integer gives a hint in the RS perturbation method. To avoid that e and \ip) 
should be obtained in pair step by step in this method, we introduce a better one where the 
effective quantum states for the system are able to share a unique set of effective operators 
such as the effective Hamiltonian and the loop current operator. 



C. Effective Hamiltonian 

The photon transition path concept leads to an easier understanding on H{e). Let us 
expand H{e) to order with the aid of Fig. [Hand Tabled To begin with, like Hs,sx 
Eq. (13T!) . each CPTP operator in H{e) holds its own identical dominant term, the order of 
which facilitates comparing its relative strength with others. According to the order analysis 
( see Appendix [A] for some details ), all of the CPTPs involving the third or higher excited 
levels, among which the one P8 Pg = -ffo,3^3,o provides the maximum correction on 0(/3^), 
can be dropped as well as the infinite weak ones bound in the three lowest levels, i.e., from 
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FIG. 4: (Color online) CPTPs of different types. The last steps of those CPTPs do not denote 
the 7io,s-like operators but the i?o,s-like ones. 

TABLE I: For typical CPTPs labeled in Fig. U their corresponding operators and the orders of 
their dominant terms on f5 are listed. The CPTPs PI, P3 and P8 denote the cases of the direct 
connection type for the first, second and third excited levels, respectively. By checking the orders, 
it is found that PI is stronger than any CPTP involving the first excited level such as P2, P5, P6 
and P7. The photon energy denominator in the self-transition operator yields that P2 is 



weaker than PI. 

Label Operator Order 

PI A = i?o,iWi,o 

P2 ^2 = ^0,lWl,l(^)^l,0 

P3 Vs = Ho,2H2,o 

P4 A = i?0,2W2,2(£)H2,0 /3' 

P5 and P6 = Vl = i?o,iWi,2H2,o 

P7 V7 = Ho,ini,2n2,ini,o p^/^ 

P8 Vs = ^0,3^3,0 



P4 to P7, and a sum of the remaining three ones PI, P2 and P3 yields one approximate 
pseudo-Hamiltonian H^^^'^^e) as 

Jrl"- ' '{e) - -Ho,0--no,lT h -Ho.l 77 -"0,2777 i4bj 

fujjLc {riWLcY ^nujLC 
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where its superscript "(3/2)" annotates that it is expanded on 0(/?^/^) and partial higher 
order terms are also included. The equation H^^/'^\e)\(fo) = ^^Iv^o) becomes a generalized 
eigen-problem 



i^f^)|y.o) = eR\v>o)+0{f3'/'), (47) 



where 



£y(3/2) _ fj Hq^iHix) iJo, 1-^^1,1-^1,0 -f^0,2-f^2,0 , .qs 

(^LC) 

Although Eq. (H71) can be solved (see Appendix O for details), an alternative but more 
general way to eliminate the e-dependence is to substitute H{e) for e in the perturbation 
terms of H{e). For instance, to deal with the term —-(j^^^Wa) , we can multiply e, a 
constant number commuting with any operator, and |<y9o) first of all, and then replace £|v?o) 
with H{e)\{po) as follows, 

- -jr T^m) = -Jir wHi^)\Vo)= - Tir w^o,oWo) + 0{/3''^), (50) 



LC) 



where only Hq^ in H{e) is kept in the final expansion. Therefore, H^^^^\e) gets rid of its 

~ (3/2) 

e-dependence but changes to a non-Hermitian effective operator H^il as 

f.(3/2) _ ^(3/2) Ho,lHl,oHo,o , . 

Generally, because this kind of substitutions can continue to increase the orders of P of the 
remaining e-dependent terms in H{e) until the result does not depend on e on the order we 
want, this approach, namely the e-H{e) substitution, can formally achieve an accurate and 

~ ~ ('3/2) 

e- independent operator Hnh which, however, loses its Hermiticity completely just like H^^j^ , 
its expansion on O^P^^"^). As discussed in the previous papers^^i^i^, the non-Hermiticity 
comes with that \ipo) is not a good effective state candidate in the equation 

Hnhl^o) = ^\^o)- (52) 

If we introduce another eigenstate {ip) with its eigenenergy being and e is rewritten as 
for the sake of symmetry, there exists an identity overlap problem as 

i^ol^o) + 5^,^, (53) 
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which is also indicated by the generahzed eigen-problem Eq. (H7|) . In fact, defining an 
operator vector Q = [Qo,Qi, ■ ■ ■]'^ (Analogously, one can also drop the ^-dependence of the 
operators Qi{e), ^2(£), • • • as we do in Appendix [B|) with its norm being 



2 



^\\=^h\J2^sSsj , (54) 

it is found that the orthogonality of the three-phase states {ip\(p) = 5^^^p in the manifold M.q 
can be expressed by the components ^nlV'o) and ^nlv^o) as 

5^,,^ = (V^ol^ii^iilv^o)- (55) 
Let us construct a new equation from Eq. fl52|) as 

HeffWeff) = e\^eff) (56) 

with the two-phase effective state \<feff) being 

l<^e//) = ^l||<^o) (57) 

and the operator H^jf being 

Heff = Q\\HnhQ^^- (58) 

According to Eq. ( l55l) . we have restored the orthogonality of the effective states. Fortunately, 
the operator H^ff is also Hermitian (see Appendix [D] for the proofs). Therefore, the effective 
Hamiltonian Hgff with \feff) can describe the manifold Aio of the three-phase system 
accurately in a compact two-phase subspace. 

We here give some comments on the availabilities of this method via a comparison to the 
UT method. According to the above definitions, it is not difficult to obtain 

Since the state \feff) is arbitrary and the vector QG^^ is explicitly unitary, our method 
exactly focuses on the manifold Aio and presents a formal solution on its corresponding 
eigenvector belonging to the transformation T in the UT method. So Hq and \(po) in the UT 
method are equivalently H^fj and Iv^e//) here, respectively. The UT method achieves the 
expanded operator S instead of T = e*'^, suggesting that it may work more efficiently when 
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the order becomes higher. The PTP approach, however, gives clear pictures to handle the 
expansions on lower orders and also successfully predicts the properties of this problem. For 
instance, since the well-known optical selection rules forbid the photons to take odd times of 
creating and annihilating processes to go back to the same level and since the corresponding 
operators a) and a are always associated with a factor proportional to (3^^^^ it is found that 
there only exist non-zero terms on the orders of (3^/^ to even powers in Heff-, some hints on 
which have been given by the orders of the dominant terms in the CPTPs shown in Table 
[H See Appendix lA] for more details. So with an arbitrary integer r, we have -^e// in the 
expansion to order (3'^/'^ as 



^Sf = Y.Heff\,^/^+0{(3'-^). (60) 

fe=0 



where the operator Heff\pk/2 is in proportion to P''^'^. In our method, the effective Hamil- 
tonian H^^ff^ may selectively keep some higher order terms for easier calculations and de- 
notations, but the non-trivial terms ife//|^fe/2 for = 0, 1, ■ ■ ■ , r are uniquely determined 
and it still bears an error on jS^"^^^^/"^ . Consequentially, one can have the eigen-problem with 
improved conditions as 

e^^'^+i) = 0, (61) 
l^r-"'^) = 0, (62) 

where k is an integer, and both e^^'' and are in proportion to j3^/'^ defined in Eqs. ( 128|) 

and fl2^ . respectively. 

Back to Eq. (IHTj) . with the method provided by Eq. and the expansion of Q\\ on 

O{0^/^) being 

(7(3/2) _ H0,lHlfi .ggV 

r(3/2) ^ f frimy 



we have the effective Hamiltonian H^j^j = yH^/j j as 



HefP = Hofi + Ho^ifi + Ho^2,o + H^^^'^\ (64) 
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where 



^0,1,0 = (65) 
= ^''^ 

£V(3/2) _ -^0,1-^1,1-^1,0 _ 1 -^0,1-^^1,0-^0,0 + -^0,0-^0,1-^1,0 /^yx 

- (/k^^c)^ 2 (/k^^c)^ ■ 

Let us scrutinize the terms in the effective Hamiltonian HfjJ^. The first term 
^0,0= (^0 1 -^tr I ^o) originates from the projection on the ground level of the oscillator. Be- 
sides including the inductance-free two-phase Hamiltonian ffo) it also takes into account the 
vacuum fluctuations of the oscillator, both of which in total read 

1 

^0,0 = -Q^C^-^Qe-J^e-^^/^^jfe^os^,, (68) 

k=\ 



where 7^ is a dimensionless factor as 



7fe = = 1,2,3 (69) 

Oik V gOLser 



with a ratio parameter 



9 = -E^ = — (70) 



Ejo ^o^cqC'o 

?co 

showing a typical Josephson energy Ejq compared to the charging energy Eco- Because the 
sinusoidal potential of each junction is equal to zero on average, the vacuum fluctuations 
equivalently flush the junction energy Ejk into a weaker one Ej/^ as 

E'j, = e-'^'^/'Ejk,k^ 1,2,3, (71) 

which indicates that the effective size of the kth junction is reduced by a factor 6"'^^/'^ which 
gives a correction maximized on 0(/9^/^). 

The non-positive term i^o,i,o providing main effects on 0{f3) relates to the interactions 
between the two lowest levels of the oscillator. The states in the manifold Aio with e <^ hwic 
hardly afford one high-frequency LC photon so that the first excited level of the oscillator 
is almost empty due to the ensuing energy punishment. Since the state occupying the 
ground level can spread into the first excited level and also accept its feedbacks due to the 
bidirectional transitions Ho,i and Hi^ between those two levels brought by the two- phase 
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flux qubit system, the almost empty excited level acts as a "mirror" for the ground one, 
which endows a correction -ffo,i,o to minimize the eigenenergy of the eigenstates. According 
to Appendix [Bl we have 

^0,1,0 = -\L{if)\ (72) 
where the /3-independent current operator 

if = C^V^sin^ (73) 

resembles Iioapif) in Eq. ([8]) with 0i, 02 and 03 being replaced by the effective phase variables 
6'i, 6*2 and ^3, respectively, and dominates in 

if = c..rj:'^^^^ne,^ (74) 

which equivalently keeps the critical current of the /cth junction modified by a fluctuation 
factor e""^*/^ like the case of E'jj^. It is worth noting that the coupling Vi0 = ^^if in 
Hint interestingly renders if and also presents the dominant terms in ^0,1 and Hifi. To 
emphasize it, we assume that those two subsystems couple with each other only by Vi0 and 
have the Hamiltonian 

iitr,4> = ilo-lL(^ify + D,ta, (75) 

where the operator 



indicates an LC-oscillator with an additional flux displacement —Lif. In a semi-classical 
picture, the above formula suggests that the average value of the current in the loop induc- 
tance is expected to be if as a function of the slow junction phases instead of a real zero 
value when L ^ 0, so if can be understood as the loop current produced by the junc- 
tions which drives the inductance to generate an additional small flux. As a result of that 
the slow-varying-function biased LC-oscillator does not change its own eigenenergy signifi- 
cantly, the inductive energy —■^L (jf^ on 0{j3) is added as one perturbation correction to 
the effective two-phase Hamiltonian, which can be explained as that the flux generated by 
if in the inductance also affects the junctions themselves. Intuitively, this self-bias effect 
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persistently lowers the potential on any point which keeps a non-zero current and always 
opposes the current direction switching. In the quantum regime, this kind of understanding 
is still supported by the facts that the estimation (y9|l)^talv^) = ^f^LC + 0{P^^'^) provides no 
effect on 0{f3) and that the inductive energy correction dominates in -ffo,i,o- Furthermore, a 
rigorous analysis in the next section also confirms that 1^°'* is the loop current operator for 
the inductance-free flux qubit. 

The term Ho,2,o shows the direct interactions between the ground and the second excited 
levels of the LC-oscillator via two-photon transitions. Photons travel forth and back via the 
bidirectional transitions ^fo,2 and H2,o, resulting in a non-positive operator 



Ejo 4y 2g ak 



T cos^A: , (77) 



according to Appendix [Bl Its main effects are on ©(Z?"^/^) contributed by the coupling 
^202 ^ in E,^^. 

Finally, the last term Ilf^'^^ corresponds to the CPTP operator ^2 = -f^o,i'^i,i(^)'^i,o 
which includes the self-transition i-ix^x[e) of the first excited level. Since the capacitive 
energy part iQg Cgp^Qg of the unperturbed Hamiltonian Hq does not commute with Hx q 
and i^o,! which turn out as functions of the effective phase variables B\ and Q^-, in the effective 
phase representation (6*1, 6*2) with 



d_d_ ■ 

deader 



=ff,forM^l,2, (78) 



simplifying the right hand side of Eq. fl67|) yields Hi^^'^^ in a symmetric form for the three 



'3 

junctions as 

(3/2) _ L \\(\Tr^-l(\ f(0)l f(0) 



ar' = I \Q^c-,;Q„i':'\ ,/r +o(,r-), (79) 



k=l k=l k=l k=l k=l 



k 



where 



/ = 41 (81) 



Therefore, the effective Hamiltonian -f^e/j taken account of four corrections of differ- 
ent types to the unperturbed one Hq. Its complicated expression indicates that treating the 
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LC-oscillator as a three-level system does not stand as an easy task on the derivations and 
analysis. First of all, unlike common perturbation situations where two subsystems couple 
with each other via a weak linear interaction, the Josephson junctions exhibiting as non- 
linear inductances keep the interaction Hint in Eq. fl26l) split into the couplings of different 
strengths. For instance, among the effective corrections in proportion to Z?^''^, Vgr donates 

one as V6(f^o|0^|^o) in -^0,0, and Vs^^ as — (^LC')~"'^V3(i7o|0^|^i)^i(^i|0|^o) in -?^o,i,o- This 
kind of terms inside the photon-transition operators Hs,si have been automatically included 
in our results while the step-by-step method should explicitly calculate them out. On the 
other hand, V20^ turns on the direct connections between the second excited level and the 
ground one, thus straightforwardly imposing the influences of this excited level without the 
help of any other excited one; otherwise, only with Vi0 its maximum feedback decreases 
to the CPTP P7 V-j = in Table [H the dominant term of which is on 

0(/3^/^). Thus our method also needs to accumulate suitable CPTPs one by one. Moreover, 
although being a part of the effective potential in iifjj^^ the operator lif^"^^ involves one 
self-transition process and performs as a correction sensitive to the eigenenergy ( see the 
pseudo-Hamiltonian H^^^'^^ie) ). This feature is not good for the analysis of the experiments 
which often alter the energy level structure of the whole system by changing the external 
flux bias. After solving the eigen-problem of H^^jJ^ ^ although the eigenvalue 5(^/2) directly 
gives 

^ = e(=^/2)+o(/52), (82) 
the effective state W^ejf) should be preprocessed as 

for further discussions. Even if the difficulties mentioned above are carefully handled, we 
should still cope with tens of terms related to ai^2,3; (3 and g. Limited by the fabrication 
conditions and other factors, the loop inductance cannot be enlarged too much, and thus 
the 0(/5^/^)-effects appear essential in rare cases. Therefore, as a compromise between 
simplicity and accuracy, we choose one effective Hamiltonian on 0{(3) rougher but optimal 
in this trade-off as 

fV(l) _ ir Hp^lHlfi , . 

H^ff - Ho,o - (84) 
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which bears an error on 0(/?^/^). Dropping the terms in proportion to or higher orders 
of f3 in Eq. (18^ yields an effective potential 



(1) 



eff 



k=l 



(85) 



identical to the one presented by Ref]l4 The corresponding normalized effective eigenstate 
approximates Iv^o) on 0{(3) as 



'eff 



with the eigenenergy e^^'^ being 



£ = £(i) + 0(/33/2). 



(86) 



(87) 



D. Arbitrary effective operator 

As mentioned above, the photon transition path method presents not only an accurate 
prediction on the eigenenergy by the effective Hamiltonian but also a full description of 
how an effective two-phase system is mapped to the three-phase one. Take an arbitrary 
three-phase operator F as an example. Assume that two arbitrary eigenstates \(p) and {ip) 
in the manifold TVIq go with their eigenenergy and e^, respectively, where \ip) may equate 
to lip). The expansions 

F = '^^\^rn)Fm,n{^n\, (88) 

m,n 

s 
s 

yield 

{^\F\^) = (^ol (j2^^rnFmA] |^o) (91) 
\m,n / 

= ii'efflFeffl'^eff), 

where we define the effective operator for F as 

Feff = J2{GmG\\'y Fm,nGnG\\'. (92) 
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The two-phase effective operator Ff.ff depends not only on -F^.n; for example, which may 
obey the optical selection rules, but also on Qrn, Qn and Q^^^ which portray all of the projected 
components of the eigenstates in the three-phase system. Especially, for F as the three-phase 
identity operator J, one can have 



[QtiGw GnG\\ ^ = hp-, (93) 



where is the identity operator for the effective two-phase subspace; for F = Htr, with 
Eq. (159!) we have a self-consistent result as 

[Ht.) = {QQ^') H,Mq' = Heff. (94) 



V. EFFECTIVE CURRENT OPERATOR 

In the three-phase system, the current operator can be achieved in different ways: on the 
one hand, the definition of the loop inductance L yields that 



$L $0 



(95) 



L 2nL^' 

on the other hand, according to Kirchhoff 's current law, the series current flowing towards 
the fcth junction(for A; = 1, 2 and 3) can be expressed as a sum of its Josephson supercurrent 
Jcfe sin 0fc and the one through the capacitor C^, which due to the time variation of the 
charge Qk is provided by the Heisenberg equation 

Qk = Htr,iQk/h = i<j>- Icksai(j)k. (96) 

One can simply find that the series current operator of the loop possesses a unique form J^. 

The virtual work principle, besides the direct derivation above, also suggests some rea- 
sonable forms as 

where Vtr is the potential term of Htr, because for the eigenstate \ip) we have 

{^|/».l^) = M||^b) = |^. (98) 

However, since the translations such as the one in Eq. (fT6|) can alter the dependence of Vtr 
on $j5s:, lots of current operators such as I^, the DC current operators IckS^^(pk for k=l,2 
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and 3, etc, are possible candidates in this approach, which all provide the same diagonal 
matrix elements {ip\Idi^^\ip) equal to de^/d^^x- Unfortunately, this method cannot inform 
us which one is proper for the non-diagonal elements. 

With the current operator being ready for the three-phase system, one can expand it in 
the oscillator-subsystem as 



{Q.m\'i^\VLn) = J3 {VmSm,n+l + V^^m+l.n) ■ (99) 

With the effective theory shown in Eqs. (!9T|) and (!92l) . it yields that 



W^l^) = -^-^J^-^ii^i^l^k^i) + i^k-i\Vk)) (100) 

k>l 



and 



$^v^(^i7^) {glGk^i + GlA) (101) 



'^^^ k>l 

where the tilde symbol labels the effective operators. 

Let us check the dependence of on the reduced inductance /?. The dimensional factor 
\j icTi, L^ ' belonging to in proportion to , indicates that this current operator 

generally diverges with the loop size. For example, the state \ip) with |((y9o|v^i)| = \ carries 
an infinite loop-current {ip\I(f\ip) when /3 — ^ 0. Oppositely, the optical selection rules zero 
out any rule-breaking term regardless of the order of its scale factor on /3, so a real dark 
state Iv^darfc) in the manifold M-^q"* where the excited levels of the LC-oscillator are entirely 
empty is forbidden to possess a current circulating in the loop due to (fiol^li^o) = 0. As 
a result, with the dimensional factor the largest term among the small perturbations 
left in the sum in Eq. ( llOip provides an inductance- independent operator which has 
been written in Eq.( |73l) and confirms that the junctions determine the loop current when 
the inductance is small enough. One interesting thing is that the effective counterpart for 
the photon number operator a^a. 



(«'«)e// = ' +om, (102) 

is mainly determined by the inductive energy L ^-^j,"^ j /2 divided by one LC-photon energy 
hjjLc- Thus the eigenstates in the manifold A^o actually look dim with the average photon 
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number being much less than one, neither dark with no photon completely nor bright with 
one or more photons. 

Expanding Eq. fllOip to order /3, we obtain the effective current operator ij^^ as 



where the LC-oscillator is involved as a three-level system and the operators in the paren- 
theses such as Qi should be expanded on 0(/5''/^) denoted by the subscript. It is clear that 
the large dimensional factor leads to deeper explorations on the ^-operators: Qi should be 
expanded at least on 0(/?^/^) in J^^^ but on 0(/?3/4) in H^^^)^ and on 0{(3^/^) in H^ff. Ac- 
cording to the optical selection rules, it is found that I^^ suffers from an error on 0(/?^/^) 
compared to I^p as 

= f;^ + 0{p'/'). (104) 

Although the direct expansion in Eq. f llOSp can be accomplished with the help of Ap- 
pendix [HI since our effective theory can cope with arbitrary three-phase operators, one can 
also achieve I^^ via applying the theory to its another definition on 0{(3). For the sake 
of clarity, Ico is utihzed as the current unit, h = 1 and e = 1/2. We construct a current 
operator 

3 

Lon = a^er^sin^fc (105) 
k=l 

where according to Eqs. (ITTl) and (fT9l) it has been divided into three components as 

3 

-^cos = a^er sin 4 COS ^^0, (106) 
k=l " 

Oik 

k=l " 



3 

4n = ttser COS 9k siu -^(j) + flpl^ (108) 

O^serP ^0Z^ g^3 + " 
k=l " 
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With the aid of Eq. (p6|) , it is not difficult to achieve that 

Q0 = Htr,iQ^/h =i^-icon, (109) 

which reads 

= (1 + fjp) - Jeos - /sin- (110) 

Applying (1 — fip) ■ to both sides of the above equation, we have 

= (1 - (jcos + /sin + Q^) + fili,p. (Ill) 

Utilizing the definition in Eq. fl92p . one can find that the effective operators of the three- 
phase ones (1 — r)^) /gm, ^IpQct, = Vi3 {I4, — honj and fj'^ict, equate to zero on 0(/?), and, thus, 
it follows that 

where 

'* — (^o|/cos|^o) = (^o|/con.|^o); (H^) 

Qf = ^ } 1 + 0{(3'^'). (114) 



[hxiLc) 

Three kinds of effects are taken into account in the above formula for The first term 1^^^ 
consistent with its definition in Eq. ( 17il) shows that the projections (i7o| cos {^serOik''^^^ 
impose the vacuum fluctuation factors e~^'il'^ to the corresponding terms in 1^^^ for /c = 1,2 
and 3. The second term —fjfBl^ is traced back to —fjpl^ in Icon as the result of the linear 
approximations afc^^a^er^ for sin ^afc^^aser^^ And the final term Q^^^ represents a tiny 
current Qcf, fiowing through the series capacitance CseT-, which can also be obtained from 
the direct expansion in Eq. (11031) or the solution presented by Appendix [El As a two-fold 
commutator, it correspondingly involves a self-transition process i-i\i{e)i-ii,Q occurring in 
the first excited level of the oscillator and, thus, explains why the self-transition processes 
are able to challenge the Hermiticity of the effective Hamiltonian. With the effective states 
\ip^eff) I'^e//)' matrix element {ipfj^\Q^J^\(pfj^) can be calculated numerically as 




{^S'/IQ;"!^'/) = bS'/) + 0(B^'% (115) 
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FIG. 5: (Color online) Features of the potential of the three-phase system in Eq. ([9]) against the 
reduced inductance /3: (a) loop current in the potential minima, (b) barrier height of the double- 
well structure, (c) distance between two neighboring minima in the potential. Other parameters 
are ai = 02 = 1, / = 0.5 and 03 is shown in the legends. 

where e^^^ and e^^^ being involved indicates that the self-transition effects distinguish the 
corresponding eigenstates by their different eigenenergy. 

Now we have a short summary of several effective current operators in the effective 
theory. The operator I^^ defined in Eq. f l73|) as an effective current operator excluding any 
inductive effect acts as the loop current operator for the inductance-free flux qubit system. 
The operator I^^ defined in Eq. ( 17^ appearing in the effective Hamiltonian H^^j^j contains 
the vacuum fluctuation corrections while both ij^^ and 1^^ treat the oscillator as a two-level 
system. The third one costing more, can include the effects brought by the second 
excited level of the oscillator and, especially, possesses a term coming from a self-transition 

~ n\ ~ (3/2) 

process which the effective Hamiltonians H^j^ and -f^e// have. Formally, the effective 

current operator like -f^e// can be truncated on arbitrary orders. However, for the next 
step, to achieve the fourth one jj,^^ accurate on 0(/5'^/^), we should no longer neglect Q^^ and 
expanding Qi to order turns out as a more cumbersome task without any surprise. In 
the trade-off between simplicity and accuracy, we choose as the optimal approximation 
for J^, which is also justified by the following numerical simulations. 

VI. NUMERICAL DISCUSSION 

To begin with, we first shortly discuss the potential of the three-phase system at the 
degenerate point 0x = ti" or / = 0x/2vr = 0.5 with ai = 0:2 = 1 on the coordinates 

(</>l, 02,03)- 

Figure [5t^a) shows the magnitude of the loop current Iq in the potential minima as a 
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function of (3, where different curves correspond to different 03. Wlien 0:3 < 0.5 and (3 
is small enough, the bias flux cannot drive a persistent current in the loop, and a zero- 
current point Pq = (0, 0, vr) mod 2tt achieves the potential minimum —ai — ^2 + c^s- As 
j3 increases, the loop begins to support a non-zero current. If is smaller than a critical 
value a3 = V2/2, which is obtained by Ig\f3=o = as-^co, the increase of P enlarges the loop 
current to achieve the value Ics = Cislco, the maximum current the loop can afford. In other 
conditions, a larger inductance always suppresses the loop current. As (3 is large enough, 
e.g. P = 5, the inductance aggressively erases the differences caused by ^3 and damps Ig 
into zero quickly, indicating the domination of the inductance f3 in this regime. A simple 
calculation shows that the phase of every junction tends to be mod 2tt and, therefore, 
the TT-phase with which we bias the circuit comes to drop on the loop inductance itself; 
for example, when (3 = 5 and = 0.8, the inductance phase (f) reaches O.GOvr. Actually, 
when f3 is large enough, a tiny but non-zero current approximately on 0{^q/2L) can make 
the inductance possess almost the external vr-phase bias and force the circuit to approach a 
possible global potential minimum — Ylk=i ^k- 

When the circuit begins to support a finite non-zero current in the potential minima, 
its direction degeneracy yields that the potential minima such as P+ = {(f)i , (f)2 , (p^) , and 
P_ = {—(pt, — 0^,27r — (p^), where (f)f23 belong to the interval of (0,7?), depart from each 
other in pair while the zero-current point Pq pins in the phase space as a saddle point of 
the barrier, thus forming a double-well potential structure. In this condition. Fig. ^h) and 
(c) demonstrate the barrier height defined by the potential difference between P+ and Pq 
and the distance between P+ and P_, respectively. It is found that a larger (3 enhances 
the potential barrier and separates further the well bottoms. Those numerical data verify 
an intuition that such a non-negligible inductance suppresses the speed of switching the 
directions of the loop current. Consequently, in the quantum regime, these properties also 
correspondingly weaken the interactions between the two persistent-current states. See Ref. 
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for a detailed discussion on the three-phase system as well as its numerical method we 
utilized.— 

To study the quantum behaviors of the flux qubit system, the tight-binding model can be 
utilized in the first step,-'^'^ and the Hamiltonian of the three-phase flux qubit with proper 
parameters can be expanded approximately in its two flux states locating in two neighboring 
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10-^ 

-4 -3 -2 -1 

FIG. 6: (Color online) Tunnel splitting of the flux qubit, in units of Ejq, vs. the reduced 
inductance f3. Parameter is selected to be equal to two typical values 0.6 and 0.8 while others 
are ai = 1, 02 = 1, g = 80 and / = 0.5. The inset with the same symbols re-scales the range of /3 
and draws the percent changes of the numerical results to the corresponding values on /3 = 0. 

potential minima respectively as 

H = 5^xlp<yz - Aa,, (116) 

where and are Pauli matrices, 5$x is the flux deviation from the degeneracy point 
$x = 0.5$05 A > is the tunneling energy between two flux states and Ip is the magnitude 
of the characteristic current possessed by the flux states. Deflne the matrix element imn 
of the current operator as imn = {iTL\I^\n) for the mth and nth eigenstates \m) and \n) 
with eigenvalues Em and £„, respectively. At / = 0.5, it yields 2A = ei — Eq and Ip can 
be calculated as the magnitude of zio- For comparison, the effective theory also gives the 
corresponding results such as 2A = ii - Eq, imn, imn-, Ip^^ = |4?l/=o.5 and Ip'^^ = |4f |/=o.5, 

~(2) 

which we utilize the tilde sign to symbolize in this section. The symbols imn represent the 
matrix elements of the effective but non-optimal current operator I^^ and imn are achieved 
by the optimal effective current operator The numerical comparisons are given as 

follows. 

Figure El plots 2A and 2A as functions of f3 based on = 0.6 and 0.8. When (3 is small 
enough, e.g., (3 < 0.1, 2A only deviates slightly from its inductance-free value. Furthermore, 
when j3 is re-scaled (see the inset), it is found that A does not decrease monotonically but 
reaches its maximum value at /3 ~ 10^'^. As mentioned above in the effective theory, the 
vacuum fluctuations on 0(/9^/^) brought by the LC-oscillator actually reduce the effective 
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FIG. 7: (Color online) Current factors Ip, ij^^ and /p"^^ in units of Ico as functions of the reduced 
inductance /?. The inset plots the deviations of the currents from the corresponding values Ip\f3=o 
on different 03 in per thousand. The parameter 03 is also chosen equal to 0.6 and 0.8 while the 
others are ai = 1, 02 = 1, g = 80 and / = 0.5. 

-3.5 

-4 
-4.5 

-5 

m -5.5 

o 

dT -6 
-6.5 

-7 
-7.5 

-8 

-4 -3.5 -3 -2.5 -2 

'o9ioP 

FIG. 8: (Color online) Relative Errors l^,"^^ and — — as functions of the reduced inductance 

p. 

sizes of the Josephson junctions, thus suppressing the barriers and enhancing the interactions 
between these two flux states. On the other hand, the self-biased inductive effects hke 
—LP/2 on 0{(3) increase the barriers and slow down the current direction switching speed. 
As a numerical order prediction, we have those two characteristic factors equal as 7^ ~ /3 
and get a critical value (3 ~ 10~^ agreeing with the data of the inset. As (3 becomes larger, 
a clear tunnel rate damping means that the self-biased effects grow up to a non-negligible 
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FIG. 9: (Color online) Magnitudes of the current matrix elements imn for the three lowest levels, 
in units of Ico^ vs. the reduced bias flux /. Two cases (3=0 and (3=0.5 are selected, and other 
parameters are ai 2 = l,a3 = 0.8 and g = 80. 




FIG. 10: (Color online) Normalized \imn\, \irnn\ and |imri| as functions of the reduced inductance 
P when / = 0.499. The data \imn\, \irnn\ and have been divided by the corresponding values 
\imn\p=o^ and other parameters are 01^2 = IjOa = 0.8 and g = 80. The inset demonstrates a typical 
error analysis. 



level. When /? > 1, 2 A is more than one order of magnitude smaller than its inductance-free 
value, and the effective result 2A decays more excessively than 2A does; in this situation, 
a small A means that two flux states of the flux qubit interact with each other weakly and 
slowly, rendering that the whole system fails to act as a useful qubit in a larger 03 such as 
as = 0.8, but as = 0.6 only makes the flux qubit slow down which may benefit the design on 
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it with a large loop inductance. It is a pleasure that when the effective Hamiltonian on 0(/3) 
fails to calculate the inductive effects that involve higher excited levels of the oscillator, the 
three-phase system with a set of traditional design parameters may no longer perform as a 
good qubit. 

To show the performances of the optimal effective current operator Fig. [7] depicts the 
numerical data of Ip^\ 1^'' and Ip vs. (3 on the cases of as = 0.6 and = 0.8. There is no 
doubt that ij^^ perfectly achieves the results in a high precision regardless of as; even when 
the inductance has a non-negligible size {(3 ~ 1), it can also correctly predict the profiles of 
the Ip — (3 curves. These curves resemble their classical counterparts Iq in Fig. [5](a), which 
infers us that the shifts of the classical potential minima introduced by a large inductance 
also take significant roles in the quantum regime. Compared to Ip and ip'^\ Ip^^ without full 
0(/5) corrections fails to describe Ip when the influences imposed by the inductance become 
notable, e.g. /? > 0.01, which also emphasizes that the effects dominate in this region. 
In fact, the inductive energy term on 0{j3) in H^ff tends to make itself minimized averagely 

~(2) 

in a relatively large (3 region which forces Ip to rise too pronouncedly to approximate the 
real value Ip. As mentioned above, the vacuum fluctuations of the LC oscillator bring in the 
0(/5^/^) effects and, thus, reduce the effective sizes of the junctions. Therefore, the currents 
are expected to decline when P is small enough to make the effects negligible, which 

is also confirmed by the inset of Fig. [71 When as = 0.8, since the net 0{(3) effects also 
depress the currents ( see Ip when j3 > 0.01 ), Ip and Jp^'' both monotonously decrease in 

—(2) 

the whole region. On the other hand, lacking full effects, Ip |a3=o.8 increases in the 

large j3 region, so there exists a minimum at /5 ~ 0.01 in the corresponding curve when the 
0{l3^^'^) and 0(/?) effects strike a balance. For as = 0.6, minima are also found to show the 
balances between the opposite 0{P^^'^) and 0(/5) effects. Both of those two types of minima 
support our previous conclusion that P ^ 10"^ ~ 10~^ is the watershed to distinguish the 
region dominated by the vacuum fluctuations. Figure M demonstrates the /3-dependence 
of the errors which A and Jp^^ bear. The linear fitting indicating that these errors are 
approximately on 0{f3^'^) sufficiently verifies our analytic conclusions. 

When a small magnitude of time-dependent flux is applied to drive the circuit, the matrix 
elements imn contribute to the strengths of the photon-assisted transition rates, significant 
for the control of the circuit. We consider its three lowest levels and plot the magnitudes 
of the matrix elements ^01,02,12 vs. the reduced flux bias / in Fig. [9l When /3 increases 
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to (3 = 0.5, the \imn\-f curves keep the same shapes approximately except for that their 
hne- widths are much narrower, meaning that it is more difficuh to control the circuit. The 
magnitudes shown in Fig. (TDl as / deviates slightly from / = 0.5, are consistent with our 
previous conclusion on the weakened interactions, and also indicate that the effective current 
operator can accurately predict the results even when (3 reaches one. The inset of Fig. [TU] 
also supports our order analysis. 

To sum up, the optimal effective current operator describes the loop inductive effects 
in good agreement with the three-phase full quantum predictions even when the size of 
the inductance L is comparable with the effective ones' of the junctions (/3 ~ 1) and, 
consequently, the circuit may perform as a less useful qubit. 



VII. CONCLUSION 



In summary, an optimal effective two-phase current operator for 3jj flux qubit has been 
obtained if considering the inductance of the circuit loop. In our classical analysis, we have 
utilized a source transformation to achieve the current form for the inductance-free two- 
phase system. Then after constructing the Hamiltonian for the three-phase system in the 
original phase space {(f)i,(j)2,4'3), "we choose a reasonable linear transformation to reformat 
the Hamiltonian in new variables, where the small inductance phase is separated as a 
single coordinate, and find that the system can be treated as the inductance-free flux qubit 
interacting with a high frequency LC-oscillator. Under the condition that the energy of 
the slow two-phase flux qubit system is small enough comparing to the LC-oscillator's, an 
effective theory has been developed for physical variable operators from the photon transition 
path method based on the BW expansion, which is also suitable for other superconducting 
circuit types. As an application for the relatively simple results which are still of high 
accuracy, the effective Hamiltonian is expanded to order P and only give an error on 0(/5'^/^). 
Besides the direct expansion, enlightened by the classical view on the circuit, we have also 
presented another simpler method in the effective theory to achieve the explicit form of 
the optimal current operator I^^^ whose corresponding error is merely on the order of 
Finally, we have verified that the optimal effective operators perfectly describe the numerical 
properties of the three-phase system. 
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APPENDIX A: ORDER ANALYSIS 

For Hint in Eq. (126|) . we define its expansion as 

= 5^VS/3'=/^ (Al) 

k>l 

where Vm]n are independent of (3. According to the optical selection rules, we have a 
(instead of j3'^/^) series expansion on Hg^s^ as 

A;>0 

where the operator vi°i is an alias for the Hamiltonian and "~ 0(/5'''~''^'''*^)" denotes that 
the dominant terms in the corresponding operator Hg ,,^ are on 

(9(^k-^i|/4)_ r|nj^g operator 

Hs,s^ combines the effects in proportion to (3^^'^^^^'^, ^|s-si|/4+i ^^^^ Gener- 

ally, the product 11^=1 ^3^,8^-1 ~ 0(/?'^^fc=i l^fc-'^fc-iD/^) can also be expanded in a 
where n and Sk are integers. 

Operators 7tIs,si introduced do not complicate the order analysis on both H{e) and Qsi^)- 
The denominator shcuLc is a scale factor in proportion to and e in Tig^g^e) is dominated 

by e'-"^ independent of (3. Therefore, one can obtain a typical term Htyp(e) in H{e) as 

where n, Sq, . . . , s„, io, ■ ■ ■ and in are integers, and Sk 7^ s^+i for k as an integer. The e-H{e) 
substitution can change H{e) into a e-independent one on a specific order via using the first 
few largest e-less terms like Hq to replace e. Mathematically, with n being an integer and 



36 



r = n/2, we can respectively expand the non- Hermit ian Hamiltonian iJ^^ to order and 
the operator Qs^^ to order as 

H^:I = H^:Iy + H%.,. +H% + --- + + OiP^^h (A4) 

= (A5) 

where Hi^i\pk/2 oc /^'^/^ and gP\ps/4+k/2 oc /J^/^+^/a for the integer k. With Eqs. dAl) and 
(1A5I) . it can be verified correct that the operator Q\\ in Eq. ( |5^ does not hold any term in 

2fc4-l ~ _____ 

proportion to (3^^ for as an integer. Therefore, the statement on H^fj in Eq. flBU]) can be 
justified without any doubt. Since the effective states and their eigenenergy are determined 
by the effective Hamiltonian Heff, we also have ^(^fc+i) ^ g g^^^ I^'q'''^^^) = for /c as an 
integer in Eqs. fl28|l and fl29l) . respectively. 



APPENDIX B: DETAILS OF OPERATOR EXPANSIONS 



In the following sections, Ico is utilized as the current unit, ^ = 1 and e = 1/2. First of 
all, the operator Hs,si ( here s may equate to Si ) in Eq. fl3T]) is 

"3 



afc cos 6k - «fc cos ( 6k + ) 



(Bl) 



fc=i 



Ho6s,sr + ttfc cos - ^ I «fce ^ ^ 



min(s,si) 



si-2t 



k=l 



t=0 



t!(s-t)!(si -t)! 



+ c.c. 



where 7^ = d -^^—^ k = 1,2, 3, have been defined in Eq. fl69|) : especially, we have 



^ V V 'k 

Ho,i = 2^afc7fce"~sin6', 



fc=i 



V2 
2 



0,2 = — X] "fc7fce ^ cos 6^/^ 
fc=i 



(B2) 
(B3) 



Equation (]B1|) as an explicit expression is consistent with the ones in our former order 
analysis such as Eq. flA2p . For s = Si, the operator Hgs holds its dominant term Hq 
independent of /?. Since e~~ = 1 + 0(/9^/^) and the largest term among (27^)'^''''^^"^* in the 
sum is on 0(/?'*~''^'/^) when t equates to min(s,si), we have the dominant term of Hs,si on 
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0(/?l* g ^ g_^_ ii ig (j^g lY^Q optical selection rules that the fluctuation factors 

9 

can be expanded in a /3^/^-series. Thus, the 



e 2 as well as the sums about (^7^ 



operator Hg^si is capable to be expanded in the same way. Equation (]B1|) also yields 

Ho,o = Hs,s + 0{(3^^^), (B4) 
where s 7^ 0, and the terms in proportion to miss due to the optical selection rules. 



We expand Gi,2,3{^) as follows, 



(B5) 



-f^i.o _^ {Hi,i — e)Hifi ^ (-^1,1 — e^Hifi ^ ^1,2-^2,0 _^ (9(^9/4^ 



LC) 



2{nujLcy 



Ho, Hi n 



Q2{e) 



2{huJLC 

n2,o + 'H2,2{e)'H2fi + 71^2,17^1,0 + 0{P^'^) 

H2,iHio 



H2,0 _^ (-?^2,2 — £)H2fi 



2hWLC 

H2.0 



+ 



Ho, H2,o 



H2,iHio 



2huJLc ^{fl^LcY 2{f}hJLc 
Hao 



where since from the effective Hamiltonian we know that 

H{e) = Ho + 0(/3i/2) = Ho,o + 
we utihze the e-H{e) substitutions on speciflc orders of P as follows: 



g-^1,0 
(#1,1 - efH,, 



{huj 



LC) 



■|<^o) 



#0) 



{Hl^ + e^-2H,^,e)H,,^ 



Hl^Hifl + HiflH{ef - 2#i,ii^i,oiJ(£) 



[huj 



LC) 



Wo) 



Ho 



Ho, Hi o 



^H2.0 



LC) 



{nujLcf 

H2flH{e 



LC 



|^o)+0(/3^/^), 

.|(^o) + 0(/5«/^). 



H2,oHo 



A{huj 



LC 



(B6) 



(B7) 
(B8) 

(B9) 
(BIO) 



(Bll) 
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APPENDIX C: SOLUTION OF GENERALIZED EIGEN-PROBLEM 



For the generalized eigen-problem, due to the perturbations, the positive and Hermitian 
operator R suggests that Eq. Wl\ can be converted to an eigen-problem as 

^_i^(3/2)^_i ^ , {r^^S)+0{(3'), (CI) 

where the initial "+0(/5^'''^)" has been improved to "+0(/9^)" due to the previous discussions 
on the optical selection rules. Expanding R to order in Eq. fICip yields 

£r(3/2)| (3/2)v _ ^L^{3/2)\ ^/^2x /p^N 

where an effective Hamiltonian -f^e// independent of e reads 

£7(3/2) _ fr Ho^lHifi -f^0,2-f^2,0 -f^O,l-f^l,l-f^l,0 /pg\ 

and the effective eigenstate \'^^elf) 

The effective Hamiltonian H^^/p is Hermitian since the transformation {^Rr^ ■ i?^^ j does 
not alter the Hermiticity of H^^'^\ One remarkable thing is that {(f^^jp) is naturally nor- 
malized on 0(/5^/^). Since the eigenstate is normalized as (v^lv^) = 1, we expand it to 
order and have 

1 = i^i^) = (v^iy/vsf ) + om. (C5) 



In sum, Eq. (IC2I) is consistent with Eq. (IMI) as an eigen-problem which covers the eigenstates 
in the manifold Aio for the whole three-phase system on 0(/5^/^). 

APPENDIX D: PROOFS OF HERMITICITY OF Hejf 

First, let us calculate the value of the operator D 

D = qn.,, - Hi,gl (Dl) 
39 



According to Eqs. (152|) and (l55ll . applying (-^ol ■ Iv^o) to Eq. (IDip yields 

(V^ol^lv'o) = {^olQiHnhlfo) - {^o\Hi,gf\\vo) (D2) 

= 0. 

Assuming the dominant term Z)*^*^) of D is proportional to /?'^''^ with k being an integer, we 
can expand Eq. flD2p to order as 

(4°)|D(^)|,.r)-0. (D3) 

As the projected components and \<f^o^) are arbitrary eigenstates of the unperturbed 

Hamiltonian Hq, it yields 

Z)(^) = 0, (D4) 

and, thus, 

D = 0, (D5) 

for k being arbitrary. It follows that 

g^^H^, = HlGl (D6) 

Finally, we achieve that 

Heff = GllHnhG^' (D7) 

= 

APPENDIX E: CALCULATIONS ON EFFECTIVE OPERATORS 

For the charge operator 



= (El) 



we have its effective operator on 0(/5) as 

/ (Xserg Ql - Ql 



2p 2 



'0 )-no 
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With lyj), and being defined in Sec. lIVDt we apply ('?/'| ■ to and expand it 

on 0(/3) as 



= i{e^ - e^){ilj\Q^\ip) 



— e 



f(2) rV 
^4, >^0 



(^1 



.2 



ol 



^2 f(2) , f(2) 2 _ 9^ f(2)^ 



LC 



tt2 ?(2) , f (2) fT2 9 fV f (2) /V 



[huj 



LC) 



(E3) 



where we utihze 



Therefore, we have 



^o|<^o)=^^|<^o) + 0(/5i/2), 
Ho\^o)=e^\iP,)+0{l3^'^). 



(E4) 
(E5) 



Q 



fT2 ?(2) , f (2) fj2 r)fT ?{2) rr 



(E6) 



which is the same as Eq. (11141) . 
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Scrutinizing the Kirchhoff 



equations in Ref. 



14 is of less difference from the classical discussion 



2d . No cross term in the classical regime mathematically is equivalent 



on Hamiltonian in Ref. 

to the corresponding quantum case due to the same matrix A making both ^AC^-^^A-^^ 
(in the classical regime) and AC^^^^^A-^ (in the quantum regime) take [0, 0, 1]"^ as one of their 
eigenvectors. 

3^ As we know, in the new coordinates (</>, 6i, 62), the three-phase potential is of a crystal-slab style 
where the periodic variables 9i and 62 represent the in-plane directions while the phase the 
normal one. Therefore, the solutions on 9i and 62 should employ the Bloch theory widely used in 
the solid-state physics. Since no real crystal boundary condition such as the Born- Von Karman 
boundary condition exists in the superconducting phase space, in a full quantum description, 
we need not consider its energy band structures in the reciprocal lattice space for most of cases 
but the size of the periodic cell in the superconducting phase space should be minimized as we 
do in order to get the eigenenergy at the original point of the first Brillouin zone. However, Ref. 



22I selects a larger cell and takes account of the tiny energy splittings according to the inter-cell 



interactions, which vanish in our numerical results. 
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